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ABSTRACT 

A nonintrusive technique for measuring dynamic gas 
density properties is described. Molecular Rayleigh 
scattering is used to measure the time-history of gas 
density simultaneously at eight spatial locations at a 
50 kHz sampling rate. The data are analyzed using the 
Welch method of modified periodograms to reduce 
measurement uncertainty. Cross-correlations, power 
spectral density functions, cross-spectral density functions, 
and coherence functions may be obtained from the data. 
The technique is demonstrated using low speed co-flowing 
jets with a heated inner jet. 

INTRODUCTION 

Quantitative techniques for measurement of fluid flow- 
properties in aerospace test facilities are needed to obtain 
performance data, to assess the performance of new 
computer codes, and to establish facility flow conditions. 
In many test programs, conventional probe measurements 
are not adequate because they perturb the flow under 
study to an intolerable degree. This is a particularly 
severe problem in supersonic flows, where physical 
probes generate shocks and cause large changes in the 
flow patterns. In turbulent flows, the presence of a probe 
will alter the turbulence characteristics. 

In this paper, we describe a flow diagnostic technique 
based on the molecular Rayleigh scattering of laser light. 
The objective of the work is to develop an unseeded, 
nonintrusive dynamic measurement technique for the study 
of turbulent flows in NASA test facilities. This technique 
will provide aerothermodynamic data that are not 
presently available. These data will be important in 
studies such as the growth and decay of turbulent 
fluctuations. The effort is part of the non-intrusive 
instrumentation development program supporting 
supersonic and hypersonic propulsion research at the 
NASA Glenn Research Center. In particular, this work is 
directed to the measurement of fluctuations in flow 
velocity, density, and temperature. 


One of the main objectives in jet noise research is to 
identify noise sources in the jet and to determine their 
contribution to noise generation. In particular, 
researchers have focused on the correlation of fluctuations 
in flow parameters with far field noise 1 1]. In this work, 
we describe a nonintrusive technique for measuring gas 
density fluctuations. Density time-history is measured 
and used to obtain the power spectrum of the density 
fluctuations. Simultaneous multiple point measurements 
are used to obtain density cross-correlations and cross- 
spectra. The technique is based on the simplest 
application of Rayleigh scattering, where only the total 
intensity of the Rayleigh scattered signal is measured. 
The gas density is proportional to the total scattered light. 
Because the Rayleigh scattering cross sections is a function 
of the molecular species that make up the flow, knowledge 
of the gas composition is generally required; however, this 
is not a concern in this work, which is directed toward air 
flows where the composition is well defined. Although 
this basic technique has been used by various researchers, 
including us, in this paper we extend the technique to the 
simultaneous measurement of density at a number of 
locations with a high frequency response (10kHz). The 
resulting measurements allow the computation of the 
density power spectrum at each of the measurement 
locations. In addition, the cross correlations and cross 
power spectra can also be obtained. An important 
practical feature of Rayleigh scattering is that because it is 
a molecular scattering process, no seeding of the flow is 
needed. In fact, the elimination of particles in the flow is 
a requirement for accurate density measurements. 

We have previously measured time-average and phase- 
average density in a screeching jet [2] and used spectrally 
resolved Rayleigh scattering to measure time-averaged 
density, velocity, and temperature in a underexpanded 
supersonic, shock-containing free jet [31. We also 
measured density fluctuation spectra in fully expanded 
supersonic jets with Mach numbers 1 .4 and 1 .8 [4]. 

In other work, LDV has been used to measure velocity 
power spectrum [5]. Rayleigh scattering has been used 
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(with LDV) to measure mean density in supersonic, high- 
temperature flow [6], Two pulsed lasers were used to 
measure spatial and temporal correlations in an 
underexpanded supersonic free jet [7j. 

We have previously examined a scheme for dynamic 
measurement of velocity and temperature using spectrally 
resolved Rayleigh scattering [8] at a single point. Our 
eventual goal is to measure space-time correlations of 
quantities of interest in basic studies of jet noise. In 
particular, we plan to measure two-point correlations of 
pV 2 (the product of density and velocity squared). These 
data will be used to relate properties of turbulent flow to 
emitted acoustic noise. 


Since Rayleigh scattering is a relatively weak process, the 
uncertainty in the measurements often is set by the photon 
statistical noise (shot noise), which determines the lower 
bound on measurement uncertainty. For example, the 
variance in the number of photoelectron counts for a 
Poisson process is equal to the mean number of counts. 
Thus the lower bound for the relative uncertainty in the 
measurement of gas density p, is equal to the square root 
of the variance divided by the mean counts. 


g(p) 

p 


U,A/ 


( 2 ) 


where At is the measurement time interval. 


In this paper, we discuss the application of molecular 
Rayleigh scattering for simultaneous, multiple-point 
dynamic gas density measurements. The multiple channel 
time history records are used to calculate probability density 
functions (pdfs), power spectral density functions, cross- 
correlation functions, cross-spectral density functions, and 
coherence functions. The effect of photon shot noise, stray 
laser light scattering, and particle scattering are discussed. 
Experimental data are presented for measurements taken in 
a hot, low-speed free air jet to illustrate the technique. 

THEORY 


Rayleigh Scattering 


For a gas with density p illuminated by a laser beam with 
wavelength A,, and power P (t , the expected count rate of 
detected photoelectrons resulting from light scattered into 
solid angle Q is 


/L 


£P„pL x A„ Q 


mhc 


da 

da 


sin- j 


0 ) 


where L x is the length along the beam of the scattering 
volume, dcr/df2 is the Rayleigh differential scattering cross 
section [9], and x * s the angle between the electric field 
vector of the (linearly polarized) incident light and the 
direction of the scattered light. Generally experiments are 
arranged so the incident electric vector is perpendicular to 
the scattering plane (i.e., “s” type) and X - 90°; f° r this 
situation the scattering is independent of the scattering 
angle. However, if the electric vector is in the scattering 
plane (i.e., “p” type), the scattering is dependent on the 
scattering angle and approaches zero at a 90° scattering 
angle. In this equation. £ is the overall collection 
efficiency (including the detector quantum efficiency), h 
is Planck's constant, c is the velocity of light, and m is the 
molecular mass. 


As an example of the expected scattering intensity, about 
40 million photoelectrons/sec would be detected by an^/4 
optical system with a collection efficiency of 5% from a 
0.8 mm long probe volume of nitrogen at a pressure of 
1 atm and a temperature of 293 K irradiated by a 5 W. 
532 nm laser beam. This high photoelectron rate means 
that independent density measurements could be obtained 
at a high rate. If we obtain 10.000 independent 
samples/sec, the mean number of detected photoelectrons 
would be about 4000, and the corresponding uncertainty 
in each density measurement would be 2%. Note that this 
lower bound represents the situation with no additional 
noise sources. In practice, one often encounters addition 
noise from ambient background light, stray laser light, 
and scattering from particles in the flows. 

Statistical description of dynamic Rayleigh scattering 

The molecular Rayleigh scattering from a turbulent lluid 
can be considered as a doubly stochastic process because it 
includes both the random nature of the fluid density, as well 
as the Poisson statistics of the light detection process ( I Of 
Because of the finite bandwidth of any real system, the 
detection process can be modeled as a sequence of impulses 
filtered by a linear filter. Consider a total detected photon 
count rate M0 = &r(0 + A H ( t) , where A is due to 
Rayleigh scattering, and A H |/) is due to other sources. If the 
light detection system frequency response is H(f), then the 
measured one-sided power spectral density function of the 
output can be written 

= [2A, +2 J R + G XK (fj\\H(f)\ 2 (3) 

where G x (f) is the spectral density function of the fluid 

density fluctuations. (Here we assume that the mean has 
been subtracted from the signal to avoid a large component 
at f- 0.) As shown in figure 1, the measured power spectral 
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density is the desired spectrum due to density fluctuations 
added to a constant level due to shot-noise. 


Another potential source of noise is a result of Mie 
scattering from particles in the ilow. The transit time of a 
particle with velocity V crossing a probe volume with width 
iv is r = w/V. For a top-hat detector time response with 
amplitude A , the two-sided spectral density function is 


s„(f)= 


( 

At 

v 


sin nfr\ 
nfT j 


(4) 


which has a 3dB point (half power) fy h = 0.44/ r. 


Rjr) = Jimyj x(t)y(t + r )dt O) 

which is related to the covariance function 

R M= C M + V^, (8) 

It is convenient to remove the mean from the data before 
processing. In this case RJ T) = CJ r). For the special case 
of x = y, RJ t) is the autocorrelation function. 

The two-sided cross-spectral density function is defined as 
the Fourier transform of the correlation function 


In addition to the direct scattering from particles passing 
through the probe volume, a smaller signal can also be 
generated by Mie scattering from particles passing through 
the laser beam outside the direct field of view ot the 
detector. The light scattered from these particles can reach 
the detection optics by means of multiple scattering from 
surfaces in the field of view. 


Sjf) = J R xy (T)e~ i2¥T dT , - 00 < / < 00 (9) 

Note that S J-f) = S*Jf) - S yx (f) 

In presenting experimental results it is common to use the 
one-sided spectrum G X} (f) = 2S VV (/), p>0 


For many applications, the desired information is expressed 
as statistical quantities such as correlation functions and 
power spectra. For simultaneous, multiple-point 
measurements one can obtain additional quantities including 
cross-correlation functions, cross-spectra, and coherence 
functions. Consider the time history of two stationary 
random processes x(t) and y(t), representing density 
measurements at two locations. The mean, variance, and 
covariance are [ 1 1 ] 

Af, = £ [4 av = £[y] 

= e [( a ' - AV ) : \<r; = - ft , )' (5) 

=£[f- ^X-v-A'y)] 

Note that these quantities do not depend on the time history 
of the fluctuations, but only on their probability distribution 
functions. 

Dynamic quantities, such a correlation functions and 
spectra, require the actual time histories of the data. We 
can then define the covariance function 

C„ ( r) = f[{x(r) - n x }{ y(/ + r)-// v }] 

7 (^) 

= 1^7! WO - AA }{ v(/ + t) - fjjlt 

where the correlation function is 


The one-sided cross-spectrum, which is complex, can 
written in terms of its real and imaginary parts 

Gjf) = 2 J RjTk '^dT = Cjf ) - iQj ffm 


where Cjf) is the coincident spectral density function 
(cospectrum), and Q xy (f) is the quadrature spectral density 
function. These can also write in terms of magnitude and 
phase 




OJJ) = tan 


QJf) 

Cjf) 


(II) 


The magnitude squared normalized by the prcxluct ol the 
spectral density function is known as the coherence function 


r ..(/)= 


G„(./T 

c..(/)c.,(/) 


(12) 


which has a range [0,1]. The coherence function is a 
powerful tool that finds use in a variety of applications. 
Cross-spectrum analysis offers the advantage over cross- 
correlation analysis in that information may be obtained 
even if the propagation is dispersive [111. 
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Welch method of modified periodograms 

We use the Welch method of modified periodograms[ 12] 
to calculate the power spectral density (PSD) and cross 
spectrum density. In this procedure, a long data record 
sampled at rate/; for time T (total samples = N = Tf s ) is 
subdivided into a number K of smaller records (which 
may be overlapping), each of length L samples. The 
modified periodgrams of each subrecord are calculated 
using a data window, which are then averaged to obtain 
the estimate of the power spectral density. The frequency 
resolution of the resulting spectral density is thus f s /L. By 
overlapping the segments by one half of their length, a 
near maximum reduction in the variance in the PSD 
estimate is achieved: the variance in the estimated spectral 
density is reduced by a factor of 1 1/9K compared to the 
variance of a single spectral estimate. The modified 
periodiogram of the k th subrecord is given by 

where f„ = n/L for n = 0 . . . L/2, and 

(14) 

for it = 1.2 K 

and 

u=j 5 > 2 0 ) < 15 > 

L /=(> 

with W(j) being the window function. (In our work we used 
the Bartlett window.) 

The spectral estimate is the average of these modified 
period iograms 

A A . =1 

The estimates of the cross spectral density are obtained in a 
similar manner. 

EXPERIMENT 

Setup 

The flow used to demonstrate the multiple-point Rayleigh 
scattering technique consisted of a 12 mm diameter 
electrically heated free jet with a 52 mm diameter ambient 
temperature co-flow (fig. 2). The mixing of two flows of 


different temperatures creates fluctuations in the density. 
The air supplied to both flows was filtered using a three- 
stage filtration system to remove particles and oil vapor. 
Stainless steel screens were placed in both flows to reduce 
flow non-uniformity. The inner heated jet exit extended 
15 mm beyond the co-flow exit plane. A type K 
thermocouple and total pressure probe were used to 
determine the mean flow properties along the jet axis. 

The beam from a diode-pumped, solid-state, frequency- 
doubled, Nd: Vanadate laser with 5 W of output power at 
532 nm was directed along the jet axis and into the inner 
nozzle, which served as the beam dump. The diameter of 
the beam at the measurement region was about 1 mm. 

The Rayleigh scattered light was collected with a 50 mm 
focal length, //l. 2, 35-mm camera lens at a 90° scattering 
angle and focused onto a multiple anode photomultiplier 
tube. This tube consists of a 32 channel linear array: with 
each photocathode having an active area 0.8mm wide x 
7 mm high. The height of the active area was limited (to 
reduce detection of stray light) to 3 mm by a mask placed 
in front of the detector. The distance between 
photocathode active areas is 1 mm. For this work, we 
used alternate photocathodes to increase the separation 
between probe volumes. The quantum efficiency given 
by the manufacturer at 532 nm is about 8 %, and the 
anode pulse rise time is 0.2 ns. 

The positions of the eight probe volumes (set by the 
magnification of the light collection lens) was measured 
by scanning a 1 mm diameter aperture, illuminated by 
diffuse laser light, through the measurement region. The 
results of these measurements showed that the separation 
between probe volumes was 3.82 mm. 

Data Acquisition 

Only eight channels of the multiple-anode PMT were 
used in this work because of the limited number of A/D 
channels in the data acquisition system. The output of 
each PMT channel was terminated in a 10 kI2 load 
resistor. The maximum allowable PMT mean 
current/channel of 10 pA is thus converted to 0.1 V. A 
personal computer based data acquisition system was 
used to digitize the eight analog signals. The combination 
of the load resistors and the PMT capacitance results in a 
limited bandwidth. Data were taken on a single channel 
at a 400 kHz sampling rate to determined the system 
frequency response. As shown on figure 3, the 3 dB 
cutoff frequency is about 10 kHz, with a rolloff of about 
20 dB/decade. Ideally, because of the slow rolloff of this 
first order RC filter, anti-aliasing filters should be used. 
However, for this work, we lacked a sufficient number of 
filters. Thus, with the 50 kHz sampling rate that was 
used, there will be a small amount of noise aliased from 


NASA/TM — 1 999-209295 


4 


ii 



higher frequencies into the 0-25 kHz frequency band. 
The laboratory was darkened during the data acquisition 
to eliminate the detection of any ambient light. Also, the 
A/D system does not actually sample the eight channels 
simultaneously, but instead samples them sequentially. 
Thus, there is a time separation of 2.5 |is between samples 
from adjacent channels. However, this should not result 
in significant errors for the frequencies of interest in this 
study ( < 10 kHz). 

The relative response of the eight PMT channels was 
determined using the spatially uniform light source 
generated by a 200 mm diameter integrating sphere. Low 
level laser light was introduced into the sphere. Figure 4 
shows the resulting power spectrum for each of the eight 
channels. For this test, channel number 1 included a 
10 kHz low-pass anti-aliasing filter to evaluate its 
performance. The other channels did not include a filter. 
A 50 kHz sampling rate was used with 200,000 samples 
taken for each channel. The response of the eight 
channels was reasonably uniform, with the standard 
deviation from the mean being less than 8%. Note that 
the response of the channel with the anti-aliasing filter 
was reduced by a factor of two relative to the other 
channels. 

Results 

Measurements along flow axis 

Data were taken along the axis of the low-speed heated air 
jet, with the imaged region being located at two positions. 
In the first region, the ratio of the distance from the jet 
exit plane to the jet diameter (X/D) ranged from 1.4 to 3.6. 
In the second region X/D ranged from 3.1 to 5.3. At each 
of these positions, eight channels of Rayleigh scattering 
data were obtained at a 50 kHz sampling rate. The data 
acquisition time was 4 seconds, yielding 200,000 samples 
for each channel. 

Data were taken at a number of How conditions, with the 
inner jet exit plane temperature at ambient temperature, 
200°, 300°, 400°, 500°, and 600° C. Here, we only show 
data taken at ambient temperature and at 600°. The 
power spectral density estimates obtained using the Welch 
method of modified periodograms described above for the 
eight channels in the two test regions are shown in 
figures 5 and 6. For these power spectra estimates, the 
subrecord length was 2048, so the frequency resolution 
was 24.4 Hz. 

Figure 5 shows the spectral estimates for the inner jet at 
the same temperature as the outer coflow. The velocity at 
the inner jet exit plane, determined from a total pressure 
measurement, was about 8 m/sec. As, expected these 
show very little density fluctuations. The small increase 


in the spectral density at low frequencies may be caused 
by secondary light scattering from particles in the laser 
beam outside the measurement region. Or, it may have 
been the result of scattering from a small number of 
residual particles in the flow. 

Figure 6 shows the power spectral estimates of the density 
for data taken at an inner jet exit temperature of 600° C 
and a exit velocity of 19 m/sec. The spectra exhibit very 
large values at low frequencies, about two orders of 
magnitude above the spectra for unheated flow. The 
bandwidth of density fluctuations is about 2 kHz. 

Information on the flow evolution can be obtaining by 
using the simultaneous multiple channel data. The cross- 
correlation between adjacent channels (such as shown 
figure 7) gives a measure of the convective velocity ol 
turbulent structures in the flow. For channels 1 and 2 
(X/D = 3.1, 3.4), the delay of the cross-correlation peak is 
0.5 msec. Using the 3.8 mm channel spatial separation 
gives a 7.6 m/sec convective flow velocity. Figure 8 
shows the convective velocity determined from cross- 
correlation for seven locations in each of the two 
measurement regions. 

More information can be obtained by calculating the cross 
spectra. For example, the phase of the cross spectrum 
between channels 1 and 2 for the 600° C flow is shown in 
figure 9. The linear slope indicates a non-dispersive 
phenomenon, and the velocity calculated from the 
0.1807Hz slope is 7.6 m/sec, the same as found from the 
cross-correlation function. The advantage of the cross- 
spectrum analysis is that dispersive phenomena can also 
be studied. 

The other quantity of interest that can be obtained from 
the cross spectrum is the coherence function (as discussed 
above). Figure 10 shows the coherence functions 
calculated for one of the measurement regions between 
channel 1 (X/D = 1.4) and the other 7 channels (at X/D = 
1.7, 2.0, 2.4, 2.7, 3.0, 3.3, and 3.6) with 600° C flow. 
Note the appearance of density fluctuations in figure 10b 
with frequencies of about 1 kHz for the coherence 
function between X/D = 1.4 and 2.0. This peak occurs at 
lower frequencies at wider separations, which is probably 
a result of the slowing of the flow as shown in figure 8. 
The subrecord length used for these coherence functions 
was 512, resulting in a frequency resolution of 97.7 Hz. 

CONCLUDING REMARKS 

This preliminary work demonstrates the potential of 
molecular Rayleigh scattering as a technique to study the 
dynamic behavior of turbulent flows. For flows with 
stationary turbulence parameters, the relatively weak 
Rayleigh scattering can be analyzed by taking long record 
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lengths and then averaging using the Welch method of 
modified periodograms to reduce the variance of spectral 
estimates. Although the work presented here used low 
speed co- flowing jets with the inner jet heated to generate 
density fluctuations, the technique can be extended to high 
speed compressible flows which have large density 
fluctuations related to the flow turbulence. Furthermore, we 
plan to extend the Rayleigh scattering dynamic 
measurements to include velocity using a previously 
developed technique based on resolving the Rayleigh 
spectrum with a Fabry-Perot interferometer. The goal is to 
measurement correlations of the product of density and 
velocity at two or more locations. This will be a powerful 
new tool for studies of the generation of noise by 
supersonic jets. It will be of particular importance because 
it is a non-intrusive technique. Use of any physical probe 
measurements in high speed flows provides questionable 
data, at best, because the flow under study is altered by the 
probe. 
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Frequency 


Fig. 1 - Expected spectrum for Rayleigh scattering 
density measurements in turbulent flow. The dc 
component of the signal is not shown. The shot noise 
floor is due to the photon statistical noise and is 
proportional to the mean signal level. 



Fig. 3 - Single channel of data sampled at 400 kHz. This 
shows that the 3 dB cutoff frequency is about 10 kHz 
when no additional low pass filtering is used. 



Fig. 2 - Experiment to measure density fluctuations at 
eight locations along axis of heated free jet. Each probe 
volume, defined by the image of the active area of the 
PMT photocathode and the laser beam width, is 1 .53 mm 
wide by 1 mm high. The separation between probe 
volumes is 3.82 mm. 



Frequency, Hz 


Fig. 4 - Power spectrum of laser light from integrating 
sphere. Channel 1 includes 10 kHz antialiasing low pass 
filter; level is reduced by factor of two for this channel 
relative to the other channels. Channels are offset by 
powers of 10. 
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Fig. 5- Simultaneous power spectral density for cold flow 
at eight axial locations (a) X/D from 1 .4 to 3.6, (b) X/D 
from 3.1 to 5.3. Spectra are offset by powers of 10. Based 
on 200,000 samples for each channel with 50 kHz 
sampling rate. 


Fig. 6 - Simultaneous power spectral density at eight 
axial locations (a) X/D from 1 .4 to 3.6, (b) X/D from 3. 1 
to 5.3. Exit temperature is 600° C. Spectra are offset by 
powers of 10. Based on 200,000 samples for each 
channel with 50 kHz sampling rate. 
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Delay, msec 


Fig. 7 - Cross-correlation of channels 1 and 2, located at 
X/D = 3. 1 and 3.4. Corresponding velocity is 3.82 
mm/0.5 msec = 7.6 m/sec. 



Fig. 9 - Cross spectrum phase for channels 1 and 2, X/D 
3. 1 and 3.4 for jet exit temperature of 600° C. Slope for 
the frequency ranage of 0 to 2000 Hz is 0. 1 80°/Hz. 
Corresponding velocity is 7.6 m/sec (same as given by 
cross correlation shown in figure 7). 



X/D 


Fig. 8 - Velocity calculated from cross-correlation of 
adjacent channels. X/D based on mean of two channel 
positions. Jet exit temperature = 600° C. 
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Frequency, Hz 


Fig. 10 - Coherence functions between X/D = 1.4 
and X/D = 1 .7 (a), 2.0 (h), 2.4 (c), 2.7 (d), 3.0 (e), 
3.3 (0, 3.6 (g); jet exit temperature = 600° C, jet 
exit velocity = 17 m/sec. 
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